Diffusion-Limited One-Species Reactions in the Bethe Lattice 
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We study the kinetics of diffusion-limited coalescence, A + A A, and annihilation, A — > 0, in 
the Bethe lattice of coordination number z. Correlations build up over time so that the probability to 
find a particle next to another varies from (p is the particle density), initially, when the particles 
are uncorrelated, to [(2 — 2)/z\p'^, in the long-time asymptotic limit. As a result, the particle 
density decays inversely proportional to time, p • 
asymptotic constant value. 



Xjkt, but at a rate k that slowly decreases to an 



PACS numbers: 82.20.-w 02.50.-r, 05.40.-a, 05.70.Ln, 
INTRODUCTION 

Diffusion-limited reactions are dominated by fluctua- 
tions at all length scales. In lack of a coherent approach 
for their analysis the simplest reaction schemes, such as 
one-species coalescence, A + A ^ A, and annihilation, 
A + A 0, have attracted the most attention 1]. In one- 
dimensional space their kinetics can be analyzed exactly 
and is found to be anomalous — at variance with the 
mean- field (reaction- limited) result: the density of parti- 
cles decays as p ~ l/Vt- In two dimensions p ~ \nt/t, 
and only in three dimensions and above does one observe 
the mean-field result p ~ 1/t 

Here we analyze coalescence and annihilation in the 
Bethe lattice of coordination number z (Fig. [T|). Our 
analysis, which is based on a generalization of the method 
of intervals 0, H, 0]j shows that even in these infinite- 
dimensional objects fluctuations play an important role, 
slowing down the rate of reactions. The correlation-free 
mean- field limit is achieved only as 2; —^ 00. Other 
reaction schemes have been studied on the Bethe lat- 
tice including trapping @, diffusion-limited aggre- 
gation 0, [1] — a process related to data compression 
and the Ziv-Lempel algorithm fs'] — and annihilation be- 
tween immobile reactants ^] . Diffusion in the Bethe lat- 
tice is also of interest to probabilists because of its close 
analogy to diffusion in hyperbolic (negative curvature) 
space [lOj. From a practical viewpoint, our system mod- 
els exciton annihilation in extended dendrimers. These 
Cayley-tree-like macromolecules are being studied as po- 
tential nanoscalc antennae, for their light harvesting ca- 
pabilities 

The Bethe lattice of coordination number z may be 
obtained by starting from a single node (at shell £ — 0) 
that is connected to z neighbors, at ^ = 1. Each of the 
nodes in shell ^ > is connected to {z — 1) nodes in 
shell £+1 and the process continues indefinitely (Fig.[T]). 
When the construction stops, at shell £ = L, say, the 
resulting graph is known as a Cayley tree. The Bethe 



lattice is the interior of an infinite Cayley tree (i.e., it 
has no boundary) . It is easy to show, for example by in- 
duction, that a connected cluster of n nodes, in the Bethe 
lattice, has exactly n{z—2) + 2 neighbors (external nodes 
connected to the cluster) , regardless of its topology. The 
case of z = 2 degenerates to a one-dimensional chain and 
has been discussed elsewhere [H, [2| ■ In the following we 
assume that z > 3. We further assume that each site of 
the infinite Cayley tree is initially occupied by a particle 
A with probability p (or empty, with probability 1 — p). 
Particles hop to one of their z nearest neighbors, chosen 
at random, at constant rate 1. Thus, in one unit time 
all particles hop once, on average. If a particle hops onto 
a node that is already occupied, the two particles dis- 
appear immediately, in the case of annihilation, or they 
merge into a single particle (at the target site), in the 
case of coalescence. 
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FIG. 1: Bethe lattice, of coordination number z = 3. 



SIMPLER MODELS AND APPROACHES 

(a) The Smoluchowski Model 

To gain intuition into the problem, let us first look at 
the Smoluchowski reaction model: a trap resides at the 
origin of the Bethe lattice, and the A particles disappear 
as soon as they hit the trap, but otherwise do not inter- 
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act with one another. This process is far simpler than 
coalescence or annihilation, yet it lets us peek into the 
nature of correlations arising because of reactions (with 
the trap). 

Let pi{t) be the density of particles in shell £, at time 
t, then, 



Poit) = 0. 



>0, 



(1) 



The first equation reflects the fact that each site in shell 
£ > is connected to just one site in shell i — 1, and to 
(z — 1) sites in shell £+ 1, while all particles hop at unit 
rate; the second equation specifies the trapping reaction 
at the origin, £ = 0, as a boundary condition. Eqs. ([T]), 
with the initial condition p£{0) — p, admit the steady- 
state solution 



Pe^P 



1 - 



1 



z - 1 



(2) 



Thus, the density of particles next to the trap is depleted. 
Right next to the trap, at £ = 1, the density is reduced 
by a factor (z — 2)/(z — 1). The missing density decays 
exponentially away from the trap. It is possible to inter- 
pret the depleted mass, p—pi, as a probability density, by 
impossing normalization: p( = {z — 2) / (z — 1)^+^. From 
it we derive a characteristic depletion distance from the 
trap, or a reaction correlation length, ~ {£) = l/{z — 2). 
At z = 2 the correlation length is infinite: indeed in one 
dimension the depletion zone near the trap expands in- 
definitely, as ^/t, and there is no steady-state. For z > 2 
the correlation is finite, but only in the limit of z — > cxd 
does one attain the ideal case of reaction-limited kinetics, 
where the correlation vanishes. 



(b) Mean-Field 

Consider now coalescence and annihilation in the 
Bethe lattice, assuming that correlations can be ignored 
completely. The neglect of correlations is strictly justi- 
fied only in reaction-controlled processes and does not 
apply to the diffusion-controlled reactions we address in 
this paper, yet it yields valuable physical insights. We 
shall refer to this limit as mean-field. 

The density of particles, p{t), is the same in all nodes, 
and since particles hop at unit rate they meet one another 
at rate p^, regardless of the graph coordination number 
z. Thus, p{t) satisfies the rate equation 



p{0)=P, 



(3) 



where A; = 1 for coalescence, A + A ^ A, where only one 
particle is removed upon each encounter, and k = 2 for 



annihilation, A + A 0, where two particles disappear. 
Eqs. ^ have the solution 



Pit) = 



1 



1/p + kt 



(4) 



so mean-field analysis predicts an inverse-linear decay 
with time, p ^ 1/kt, in the limit of t — *■ oo. 

We see that the neglect of correlations is disastrous 
for the case of z = 2 (one-dimensional chain), where 
for diffusion-limited one-species reactions p ^ l/v^, in- 
stead of the mean-field prediction of 1/t. For z > 2 
the Bethe lattice is infinite-dimensional: the number of 
nodes within £ shells increases exponentially with £ — 
faster than any power (dimension) of £. Naively, one 
would expect the mean-field limit to apply then, but the 
analysis of the Smoluchowski model, above, suggests that 
correlations might play a certain role, so long as z < oo. 



EMPTY INTERVALS APPROXIMATION 

We now generalize the method of empty intervals, in- 
troduced for the exact solution of coalescence on the line, 
to the case of coalescence in the Bethe lattice. The gen- 
eralization is not rigorous, but yields reasonable results. 
Let En(t) denote the probability that an n-cluster con- 
tains no particles at time t. The particle density (i.e., the 
probability that a site is occupied), is 



pit)^l-E,{t) 



(5) 



At time t = the occupancy of nodes is independent 
of one another, so the probability that an n-cluster is 
empty, is 



E,,{o) = {i-py 



n= 1,2, 



(6) 



Consider now changes to En{t). Any motion of par- 
ticles within an n-cluster would not alter En, since the 
(non-empty) cluster cannot become empty in this way. 
En changes only when a particle hops into or out of the 
cluster. To follow these changes, we require the probabil- 
ity Fn for having an n-cluster that is empty, connected 
to a node that is occupied. We may then write 



at z 



(7) 



The first term on the RHS denotes the event that a par- 
ticle hops out from just inside the n-cluster, leaving the 
cluster empty, thus increasing The starting config- 
uration occurs with probability Fn-i- The second term, 
proportional to _F„, accounts for a particle just outside of 
an empty n-cluster that hops into it: the cluster is not 
empty anymore and £"„ decreases. Since the n-cluster 
has n(z — 2) + 2 neighbors, and each particle hop (from 
one specific site to another) occurs at rate 1/z, the overall 
rate for all of these processes is [n(z — 2) -I- 2]/z. 

In order to close the evolution equation ([7]), we need 
to express F„ in terms of the £"„. We note that the 
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event that an n-cluster is empty (probability En) is the 
sum of the events that a neighboring node is also empty 
(probability En+i) and the event that the site is occupied 
(probability Fn), see Fig. [5J Hence, 

Fn = En — En+1 ■ (8) 

Using this in ([7|), we finally obtain 

d ^ , ^ n(z — 2) + 2 , „ „ s 

-En{t) = ^ {En-1 - 2En + E^+l) . 9 

at z ^ ' 

This equation is valid for n = 2, 3, . . . A separate analysis 
for the case of n = 1 shows that it may be subsumed in 
Eq. ([HI), with the understanding that 

Eo{t) = 1 . (10) 

If the system contains a finite density of particles, ini- 
tially, then an infinitely large cluster can never become 
empty, thus 

lim En{t) =0. (11) 

n — >oo 

To summarize, we need to solve Eqs. (|9]), for n = 1,2,..., 
with the boundary conditions (fTU]) . pT|) . and the ini- 
tial condition The particle density is then obtained 
from ©. 



n-cluster 





(a) (b) 

FIG. 3; Probability F2, in the Bethe lattice of z = 3. The 
expression F2 — E2 — E3 is correct for the configuration (a), 
but not for (b), where the empty sites do not form a connected 
2-cluster (cf Fig.[2l). 

(Fig. [3]). The approximation amounts to assigning the 
same weight to configurations (a) and (b). Since there 
are four ways for the particle to leave the 3-cluster in 
configuration (a), and only one way in configuration (b), 
the difference between the actual probability of (b) and 
the approximation (E2 —E^) matters only one fifth of the 
time. Moreover, the analysis of the Smoluchowski model 
shows us that the correlation length is of the order of 
l/(z — 2), that is, no more than a link's length. Thus, 
the actual difference between the probabilities of (a) and 
(b) cannot be large. The error involved in processing 
larger clusters is of even lesser consequence. 
• The approximation's Eq. ^ is simple, linear and ana- 
lytically tractable. The simplest alternative would be a 
Kirkwood-like truncation of the infinite hierarchy of rate 
equations for clusters of size n. Note, however, that an 
n-cluster has 2" states, so there'd be 2" — 1 variables as- 
sociated with it, and the truncation yields a non-linear 
equation at the bottom of the hierarchy. 



FIG. 2: Probability F„, for finding an empty n-cluster (empty 
oval) connected to an occupied neighbor (•). The second 
symbol from the top denotes E„+i, the probability for finding 
an (n -I- l)-cluster empty. The sum of the two events equals 

Fn. 

The above argument is not exact because the relation 
Fn = En — En+1 IS correct only when the particle is at- 
tached to the n-cluster through a single link. The particle 
might be embedded in an (n + l)-cluster, splitting it into 
two (or more) pieces. In this case we cannot express F„ 
in terms of the En- In Fig. [3] we illustrate the problem 
for the Bcthe lattice of z = 3, for the case of 3-clusters 
(the smallest cluster size where the difficulty arises). 

We argue that, despite the inaccuracy, there are good 
reasons for pursuing the empty interval approximation: 

• For z — 2 the equations reduce to those found for the 
line, where the method is exact. It may be useful to think 
of 2; as a continuous variable (in the spirit of analytical 
continuation), in which case one expects good results for 
z = 2 + e. 

• The approach is exact in the limit z ^ 00, where cor- 
relations do not play any role. 

• The approximation is good even for intermediate val- 
ues of z. Consider, for example, the case of z = 3 



Results for Coalescence 

The empty interval equations may be tackled ana- 
lytically. For example, following a Laplace transform 
in the time variable the resulting system of difference 
equations with linear coefficients can be solved in closed 
form [T3| , then the transform might be inverted (see Ap- 
pendix A). It is easier to focus directly on the limiting 
behavior. For short times we assume the Taylor series 
expansion, 

00 

En{t) = a^rk^ . 

m=0 

Substituting in ^ and collecting like powers of t one 
obtains a recursion relation for the a^'"-' in terms of the 
^(m-i)^ Since an ^ = En{0) (to satisfy initial conditions) 
is known, successive coefficients can be generated me- 
chanically up to any desired order. (This does not pro- 
vide a complete solution, for convergence of the Taylor 
series is limited to a finite radius, tc = z/(z — 2).) The 
first few terms in the particle density are 

p(t)=p-ph + p^^-^P±P^t' + --- (12) 



4 



It is interesting to compare this result to the mean-field 
prediction, p{t) — p/{l +pt), that expands to 



Pit)=P 



Mean-Field. 



We see that the first two terms are identical, but that the 
quadratic term in (I12p is larger than that in the mean- 
field expansion. This reflects the fact that the initial 
particle distribution is Poissonian and devoid of correla- 
tions. Correlations, however, do develop, resulting in the 
slowing down manifest in the quadratic term. 

The long-time asymptotics is best analyzed by passing 
to a continuum limit. As the process evolves, the den- 
sity of particles drops and n need be very large to find 
any particles within an n-cluster. Under these condi- 
tions it makes sense to regard n as a continuous variable 
"length," X = na, where a is the length of a link. We 
then replace En{t) with E{x,t) and Eq. ([9]) becomes 



d 



E{x,t) 



z-2 92 



vx-—^E{x,t), (13) 
ot z ox'^ 

where v is the average speed for a particle to traverse 
a link. The boundary conditions are now E(0,t) — 1, 
lima;^oo E[x, t) = 0, and the particle concentration, 
c{t) = p/a, is, according to (0, 

dE, 

c{t) = --g^\x=o ■ 

To reconnect with the discrete limit, we take a — 1 and 
v = l. 

We look for a scaling solution, of the form 



Eix,t)^^j^) 



as i ^ 00 



Putting this back in we see that we must have (3=1 
and that $(x/t) satisfies the equation 



z-2 



where primes denote differentiation with respect to the 
variable ^ — x/t. The boundary conditions determine 
the solution 



$(0 = exp(-— ^C) 



(14) 



and it follows that 



.z-2 



p{t)^l/{ 1) , ast^oo. (15) 

Comparing the result to the mean-field decay rate 
of p ~ 1/t, when dp/dt = —p^, we see that the reaction 
rate is effectively diminished by a factor {z — 2)/z. Thus, 
the probability to find particles next to one another is 
{z — 2)/z times smaller than p^. This is reminiscent to 
the depletion of the steady-state particle density near the 
trap in the Smoluchowski model, and is explained by the 
reaction mechanism: particles in proximity are disfavored 



in the long run, because they react more promptly and 
are removed from the system. When z —^ 00 the correla- 
tions disappear and the mean-field limit is regained. For 
z — !■ 2, the case of a linear chain, the effective mean- field 
rate vanishes. This is in accord with the known result 
p ^ 1/Vt [H, which is equivalent to a radically slower 
reaction rate, of effective order three: dp/dt = —k'p^. 



Annihilation 

We now turn to annihilation, A + A — > 0. Instead of 
En we work with G^, the probability that an n-cluster 
contains an even number of particles 01 . One obtains 
for the G„'s the same rate equations as for the i5„'s, 
Eq. but with the boundary conditions; Go(i) = 1, 
lim„_,oo G„(t) = 1/2. The initial condition is too some- 
what different: The initial probability that an n-cluster 
contain an even number of particles is 



L«/2J 



n — 2m 



i(p+(l-p))" + i(-p + (l-p))" 



or. 



G„(0) = i + ^(l-2p)", n=l,2,... (16) 
The short time analysis yields 



p(<) ^p- 2p^t + 2p^ i-^ '-t^ + 



(17) 



compared to the mean-field result, p{t) = p/{l + 2pt), 
that expands to 



p{t) ^p- 2pH + ApH^ + 



Mean-Field. 



The early buildup of correlations shows for the first time 
in the second-order term and depends upon the initial 
density p: the correlations slow down the rate of decline, 
compared to mean-field, if p < 1/2, and accelerate the 
decline if p > 1/2. For p = 1/2 the effect of correlations 
shows up first only in the third-order term. This effect 
parallels the exact result in one dimension (z = 2). 

For the long time analysis we follow a similar procedure 
as for coalescence, this time leading to [l^ 



p{t)^l/{2'-^t) 



as t 



00 



(18) 



Again, in view of the corresponding mean-field equation, 
dp/dt — —2p^, the long-time asymptotic result may be 
understood if the likelihood of pairs of occupied sites is 
effectively reduced by the factor [z — 2)/z. 
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DISCUSSION 

In conclusion, we have analyzed the kinetics of 
diffusion-limited coalescence, A + A ^ A, and annihi- 
lation, A + A 0, in the Bethe lattice, by means of an 
approximation based on the method of empty intervals. 
We find a result that is mean-field in character, p ~ 1/kt, 
where the reaction rate undergoes a slow crossover from 
fc = 1 (for coalescence) at early times to k = (z — 2)/z 
at late times. (For annihilation these values are dou- 
bled.) The decrease in the reaction rate stems from a 
buildup of correlations that disfavors particles in close 
contact. In the case of annihilation, there is additional 
behavior: if the initial particle density is larger than 1/2 
then the reaction rate increases above the early mean- 
field rate, early on, before settling to the lower value of 
the long-time asymptotic limit. The buildup of correla- 
tions weakens as the coordination number of the graph, z, 
increases, and vanishes altogether in the limit of z — s- oo, 
where the mean-filed limit is regained. The correlations 
are strongest in the singular limit of z = 2, corresponding 
to a linear chain, where the density decay dramatically 
slows down to a different functional form, p ^ l/Vi [H. 

Despite the approximation involved, we believe that 
our results reflect the essential character of the true kinet- 
ics. The long-time behavior can be explained by the fol- 
lowing heuristic argument: A particle at shell £ from the 
origin advances to shell £ + 1 with probability (z — l)/z, 
and retreats to ^ — 1 with probability 1/z. Thus, the 
particle explores, on average, (z— l)/z— l/z = (z — 2)/z 
new sites per unit time. Assuming a uniform background 
of particles, this leads to dp/dt — — [(z — 2)/z]kp'^ {k = 1 
for coalescence and 2 for annihilation), in accord with 
the approximation results. (Note, however, that this ar- 
gument fails to capture the early time behavior predicted 
by the empty interval approximation.) 

Our method relies on translation invariance and does 
not generalize easily to finite lattices (Cayley trees). In 
a finite tree of L shells the particles experience a drift 
from the center to the rim, which consists of roughly 
half of the nodes already for z = 3 (the case relevant 
to dendrimers). Thus, the finite-size effect could be sub- 
stantial and must be reckoned with. Indeed, in computer 
simulations [H [il the density converges quickly to the 
long-time behavior predicted by our analysis, however, 
this lasts only until time t L (see, for example, Fig. 5 
in [13])' ^"^d thereafter the decay is considerably slower, 
p ^ l/t^-^ pH Theoretical analysis of the long-time 
asymptotics in Cayley trees remains an important open 
challenge. 
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APPENDIX A: SOLUTION OF THE RATE 
EQUATIONS 

We wish to solve the system of equations 

^E^it) = + ^ - 2g„ + g„+0 , (Al) 



with the boundary conditions 

Eo{t) - 1, 
lim S„(t) = 0, 

n — >-oo 

and the initial condition 

i?„(0) = (l-p)". 
First, Laplace-transform in time; 

Ms) = / e'''En(t)dt, 



(A2) 
(A3) 



(A4) 



"'0 

so that the problem becomes 

4>n+l - Pn4>n + 0n-l = C„ (A5) 

wherep„ = 2-h;j(5^f)q:2, c„ = - J^f^^l^ , and (jjQ = 1/s, 
4>oo — 0. Rewrite the underlying homogeneous equation, 

An+l ~ Pn An + = 0, aS 

{n + a)A„+i - 2{?2 + b)An + (n + a)Ai-i = , (A6) 

where a — 2/(z — 2), b — {2 + sz/2)/{z — 2). Upon making 
the substitution 



An = / x''~^f{x)dx, 



(A7) 



the terms linear in n may be reinterpreted as derivatives 
of f{x) (using integration by parts), and yielding the con- 
sistency relation, 

f'{x) 2 zs 2 



f{x) {z-2)x (z-2)(l-x)2 l-x' 
leading to 

^2/(.-2) 



{z~2){l-x)J' 



(A8) 



(up to an arbitrary constant factor). This completes the 
solution of the homogeneous equation. 

To solve the original inhomogeneous equation, make 
the substitutions 



6 



yielding 



This concludes the formal solution of the problem. 



— An-lKi + A„_|_i/l„-|_i — C„ , 

This has the solution 



B„ 



1 

AqS 
An~lA 



k=l 



Ak 



k=l 



An-lA, 



-Ck ■ 



(A9) 
(AlO) 



Note that the first term on the RHS of (|A9p guaran- 
tees fulfillment of the boundary condition (j)o = 1/s. hi 
is determined from the remaining boundary condition, 
4^00 — 0; and taking into account that lim9^_^oo 



00 A: — 1 

^ Ak-iAk' ^ Ak-iAk 

k=i 1=1 " i " fc^i 



^'^'-/E^. (All) 



Manageable expressions may be obtained for the spe- 
cial case of p = 1 (all lattice sites are occupied, ini- 
tially), where c„ = 0, n > 1, leading to hi — 0, and 
01 — Ai/sAq. Finally, upon evaluating the pertinent in- 
tegrals, we find for the Laplace transform of the density, 
p{s) = 1/s- (j)i{s): 



p(.) = -(l- 



U{c,2,cs) 



z-2U{c~ 1,2, cs 



(A12) 



where c — z/[z — 2) and U{-) is Kummer's confluent hy- 
pergeometric function 15]. Inversion of this expression in 
the limits of short and long times agrees with the results 
derived in the paper by independent methods. 
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